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Zurek has provided simple causal bounds for the onset of phase transitions 
in condensed matter, that mirror those proposed by Kibble for relativistic 
quantum field theory. In this paper we show how earlier experiments with 
annular Josephson tunnel Junctions are consistent with this scenario, and 
suggest how further experiments might confirm it. 
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1. INTRODUCTION 

As the early universe went through a period of rapid cooling it changed 
phase several times. A signal of such transitions would be the direct or in- 
direct observation of topological defects (e.g. magnetic monopoles) arising 
from the inhomogeneities of the ordered phase. Such defects appear because 
the correlation lengths £(i) of the order parameter fields are necessarily finite 
for transitions implemented in finite time. This is even though, for the con- 
tinuous transitions that we have in mind, the equilibrium correlation lengths 
^ eq would diverge at the critical temperature if there were time enough (the 
adiabatic approximation) . 

Kibble Si observed that this simple causality imposed useful constraints 
on domain growth and the density of defects at the time of their formation. 
Unfortunately because of our lack of knowledge about the details of the 
transitions in the early universe it is impossible to find reliable predictions 
for current observations. However, since notions of causality are not specific 
to the relativistic quantum field theory (QFT) appropriate to the early uni- 
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verse, Zurek B@ argued that similar causal bounds were valid in condensed 
matter systems for which direct experiments on defect formation could be 
performed. These bounds can be understood in several ways. For our pur- 
poses the most useful is the following: 

As an ordering transition begins to be implemented, from time t = 0, 
say, there is a maximum speed c(t) at which the system can become ordered. 
This leads to causal horizons, outside of which there are no correlations. 
Zurek proposed that the earliest time (the 'causal time') at which defects 
can form is when a local causal horizon is large enough to accommodate a 
single defect at that time. Since the defect size £(t) = £(T(t)) depends on 
temperature T, which depends on time, the causal time i is determined from 

f(t)«2 fdtc{t). (1) 
J o 

This time is not to be confused with the earlier time, that depends on the 
global structure of the system, when the coherence length becomes smaller 
than the system. The two successful experimentsB'H on defect production in 
superfluid 3 He give very clear evidence for the primacy of the causal horizon 
in determining when defects form. 

There are other ways to formulate causalityH, most simply by requiring 
that £(t) cannot grow faster than c(t). This can be imposed either before the 
transition, or after. In simple systems all bounds agree, up to numerical fac- 
tors approximately unitya. To the level at which we are working (better than 
an order of magnitude, but to a factor of a few), they are indistinguishable. 

Whatever the case, although the causal bounds are robust, the extent 
to which they are saturated depends on the details of the microscopic dy- 
namics. From the microscopic level, causality along the lines above is not 
explicit, although encoded in the relevant dynamical equations. The picture 
is rather one of order being established through the growth of the amplitudes 
of long-wavelength instabilities. The earliest time at which we can identify 
defects from this viewpoint is when the order parameters have achieved their 
equilibrium magnitudes. Qualitatively, for simple models this time is also 
in agreement with the causal time above, showing that the causal bounds 
are approximately saturated. There is no real surprise in this. The causal 
time and distance scales t and £ = £(t) are just as we would expect from 
dimensional analysis (in the mean-field approximation) and unstable modes 
grow exponentially, whereby the dependence of the causal time (and corre- 
sponding defect density) on the microscopic parameters is only (square-root) 
logarithmic. However, before the causal time we now have a picture in which 
there is a fractal thermal fuzz of potential defects, whose density depends on 
the scale at which we look. By the causal time some of these have matured 
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into the (scale-independent) defects that we see subsequently. This is seen 
in numerical simulationsB, and see the article by one of us (RR) in these 
proceedings^, where this is discussed at greater length. 

The experimental confirmationBEl of the Kibble-Zurek predictions for 
vortex production in superfluid 3 He, at better than order-of-magnitude level, 
reinforced the hope that, if we understood why the bounds were approxi- 
mately respected there and elsewhere, we might have a better understanding 
of whether they would be respected in QFT. While 3 He has an unparalleled 
richness in its order parameters i and the topological defects that can be 
produced in it, it is a difficult system in which to understand the funda- 
mentals of defect formation. As an independent test of his assumptions 
Zurek □ also proposed experiments to measure spontaneous superflow in 4 He 
and spontaneous flux generation in superconductors (essentially topologi- 
cal charge density) in much simpler one-dimensional annular geometries, for 
which predictions are easier to make. In practice, correlation lengths in both 
these systems are so small at the relevant timescales that it has been diffi- 
cult to make reasonably 1-dimensional systems, and the experiments have 
not been performed. 

It is with this in mind that, in this paper, we shall examine the spon- 
taneous generations of defects ('fluxons') in annular Josephson tunnel Junc- 
tions (JTJs) for which correlations lengths are sufficiently large that effective 
1-dimensional systems can be fabricated and the Zurek predictions can be 
tested. One of us (RM) has conducted such experiments Ii9 in JTJs in the 
context of a different experimental programme that was completed before 
the suggestion was made that the Zurek bounds were relevant in that case. 

Although the choice of materials was not optimised for the Zurek sce- 
nario the results were very striking. Firstly, defects were, indeed, detected 
on quenching the system. This is in contrast to their non-observance in some 
other experiments, on planar superconductors^ and@ A He, specifically de- 
signed to check the Zurek picture. Secondly, although the circumference of 
the ring was at least three orders of magnitude larger than a fluxon, their 
production was rare. A single fluxon only appeared a few percent of the 
time, two very rarely, and more than two, never. At the time there was no 
understanding as to why it was so difficult to produce fluxons and no way to 
estimate the likelihood of their appearance. For the purpose of the exper- 
iments this did notmatter, since those that were produced were then used 
for other purposest3. We aim to show that these experimental results are 
compatible with Zurek's predictions which, to date, provide their only quan- 
titative explanation. However, because of the compound nature of JTJs, 
causal bounds and their saturation is more complicated. 

The possibility that the Zurek-Kibble analysis could be applied to JTJs 
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was first considered brieflyEl by two of us (EK and RR), and later a more 
substantial analysis was reported elsewhere 0. This paper provides an ex- 
tended and updated presentation of the theory and future experiments. 



2. JOSEPHSON TUNNEL JUNCTIONS 

An annular JTJ consists of two superimposed annuli of ordinary su- 
perconductors of thickness d s , width Ar, separated by a layer of oxide of 
thickness d ox , with relative dielectric constant e r . We consider a quench in 
which the (spatially uniform) temperature T(t) of the JTJ varies in time as 

At tq v ' 

in the vicinity of T c = T(0), tq being the quenching time. 

Above the transition, at temperature T c , we have two normal conductors 
but, after the transition has been completed, the complex order parameters 
= Pie tlpl and ^ = P2^ lip2 of the individual superconductors are charac- 
terised by the phases ip\ and if 2 alone. The Josephson current density at 
any point x in the annulus is 

J = J C (T) sin y?(x), (3) 

where ip(x) = ipi(x)—ip2(x), and J C (T) depends on the nature of the junction. 
All other things being constant, the larger the resistance of the oxide, the 
weaker the coupling between the superconductors and the smaller the critical 
current J c (0). 



2.1. The individual superconductors 

Before considering the effect of a quench on a JTJ it is helpful, for 
purposes of comparison, to consider the effect of the same quench (|2|) on each 
of the individual superconductors, in the limit (vanishing critical current J c ) 
in which they are not connected. We also take the limit in which their 
thickness is smaller than any other scale. 

This has been analysed alreadyS by Zurek, and we repeat his conclu- 
sions. Consider a single superconducting ring with scalar order parameter 
ip — p e lL P _ We are being simplistic here in our evasion of the gauge properties 
of the theory^. However, for strong second-order transitions it should be 
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true - essentially we are looking for a benchmark. In a mean-field approxi- 
mation the relevant adiabatic correlation length has the behaviouri 

m = am) = (4) 

in terms of the fundamental length scale £o (the size of a cold defect, among 
other things). The maximum speed at which the field can order itself is 

c(t) = c(T(t)) = C <P-, (5) 
V T Q 

where Co = £o/to is given in terms of tq, the relaxation time for long wave- 
length modes. [We save the symbol cq for the speed of light in vacuo.] Note 
that c(t) shows critical slowing down. 

Imposing the causal bound gives a causal time i s for individual 
superconductors of, approximately, 

t, = V T o T Qi ( 6 ) 
at which time the causal horizon is 

ts = m) = toM=u^- (7) 

V t s T 

At these very early times the single length £ s serves to set both the 
coherence length in p, as it achieves its equilibrium value, and in <p. With 
each change in (p of 2ir will be associated with a unit of flux within the 
annulus. At the risk of causing some confusion, we could also term these 
fluxons. On average, the number N of such fluxons will be zero. However, 
the variance in N will be non-zero. If, as a first guess, we assume that the 
phase <p takes a random walk around the annulus, circumference C = 2irf, 
then for C » f s , 

W^-O&M). (8) 



2vr \27rV& 

There are many qualifications to (|8|), but this is a good starting point. 

Numerically, Zurek proposed a system-independent estimate for tq from 
Gorkov's equation as tq = 7rh/16kBT c ~ 0.15 ps for T c ~ 10 K. The end result 
is that, for tq in seconds and T c in degrees Kelvin, 



tq/T c M. (9) 

For tq = Is and T c = 10 K, t s is a fraction of a microsecond. Further, with 
a reasonable choice of £o ~ 100 nm, we have £ s ~ 0.1 mm. 
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When we first tried to understand whether the Zurek-Kibble bounds 
had any relevance to the experiments of Ref. 1C , two of usEl made the naive 
approximation that, up to the formation of fluxons, we could treat the su- 
perconductors as independent, with uncorrelated phase angles ipi and 922 • 
In that case, (A99) 2 = (Aipi) 2 + (Ay^) 2 - For a ring of circumference 0.5 mm, 
as in the experiments, we would predict AN = 0(1) from (@). This is an 
order of magnitude too large, but there is a certain amount of latitude in all 
these estimates, and it did not seem totally unreasonable. 

However, what made the approximation untenable was that, on keeping 
the superconductors otherwise identical, but changing the coupling between 
them, no fluxons were seen at a level at which they would have been seen, 
had the superconductors behaved independently. We therefore have to turn 
to coupled junctions, where we find a reason for this. 



2.2. Josephson tunnel junctions 

In order to establish causality bounds for JTJs we need to identify the 
velocity that establishes the size of the causal horizons and the length that 
characterises the fluxon size. 

The strength of Zurek's predictions lies in the assumption that, in gen- 
eral, these can be read off from the time and space derivatives alone in the 
equations of motion in the adiabatic regime. To see how well the bounds 
are saturated requires further knowledge, but we shall assume that they will 
be well saturated here, as exemplified by the 3 He experiments. Further, in 
this scenario the adiabatic approximation is pushed to its limits, by which 
we take it to be approximately valid from the causal time onwards. 

If we ignore dissipative effects from quasiparticle tunneling and surface 
losses in the adiabatic regime at temperature T < T c , 99 satisfies the one- 
dimensional Sine-Gordon (SG) equation^ 

d\ 1 av 1 . fm , 
^-¥tTjm^ = wF) smip > (10) 

provided the width Ar of the annulus, of radius f, satisfies Ar <C f and 
Ar <C Aj(T), the Josephson coherence length. In this case x measures the 
distance along the annulus, with periodic boundary conditions. Unlike the 
case of superfluid 4 He and annular (single) superconductors, it is not difficult 
to implement Ar <C Xj(T), as we shall see. 

As with other established models of defect formation, the classical equa- 
tions are only valid once the transition is complete. We shall not use the SG 
equation, together with its dissipative and other terms, to study the evolu- 
tion of defects. In fact, we do not study their evolution at all. However, 
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Eqs. 1C and 12 are sufficient to enable us, in the spirit of the Zurek-Kibble 
scenario, to identify Aj(T), diverging at T c , as the equilibrium correlation 
length £ eq to be constrained by causality. Further, the Swihart velocity c(T) 
(with critical slowing down at T = T c \. measures the maximum speed at 
which the order parameter can changeO'E3. We shall see later that, for- 
mally, c(t) has a behaviour similar to c(i), 

c(t) = c(T(t)) = c J—. (11) 
V T Q 

Away from the critical temperature dissipative effects are small. 

The classical topological defects (the JTJ fluxons) are the 'kinks' of the 
Sine-Gordon theory 



<p±{x,t) = 4arctanexp 



x — ut 



AjCZVI - (u/c(T))* 



(12) 



where \u\ < c(T). 

In fact, ( |T2| ) is oversimplified on many counts. For example, the max- 
imum speed of the JTJ fluxon is reduced in the presence of an external 
magnetic field (which reduces its behaviour to that of a damped pendulum). 
Nonetheless, at our level of approximation, the Swihart velocity remains a 
good candidate for determining the causal horizon, after very early times. 

At very early times after the beginning of the transition (t > 0), before 
the individual superconductors have the ability to superconduct, the causal 
bound is, most plausibly, driven dissipatively by c(t) of (H). If it were the 
case that Co and cq were equal, we would be happy to use the formula (|ll|) all 
the way down to t = 0. Although cq depends on the nature of the junction, 



rather than just the individual superconductors, for the junctions of Ref. |10|, 
Co ~ 10 7 m/s for the samples used. This is an order of magnitude larger 
than Cq, as diagnosed earlier. However, we note that the size of the causal 
horizon at the causal time only depends on the scale of the velocity to its 
quarter power. To the level at which we are working (better than an order 
of magnitude, but to a factor of a few), we can take c(t) as the velocity in 

©• 

Our horizon bound for the earliest time t that we can see classical JTJ 
fluxons is when the causal horizon is big enough to accommodate one. That 
is, 

Xj(t) rs 2 / dtcit) (13) 
Jo 

at t = ij. This is consistent in that, at time ij, classical JTJ fluxons that 
have been formed will be contracting at the Swihart velocityti-3. To have 
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attempted to produce them before (as would happen with a slower velocity) 
would make them want to contract faster than causality permits. 

We shall give details later, but the end result is that tj > i s , at our 
level of approximation. We now understand better why our earlier assump- 
tion of independent superconductors is not correct. As the superconductors 
establish themselves, thermal fluctuations lead to the production of phase 
rotations in each superconductor (as discussed in Section 2.1). However, 
with the conductors connected, there will be continuous cancellation be- 
tween them. Flux lines of opposite direction can annihilate in the oxide and 
disappear from the system to leave only the residual locally unpaired flux 
lines, now classical JTJ fluxons with correlation length Xj. They will be 
reduced in number to our earlier guess, in agreement with experiment. 



3. SYMMETRIC JTJs 

After these generalities, we need a specific model. For simplicity we take 
a symmetric JTJ, in which the electrodes are made of identical materials 
with common critical temperatures T c and assume a quench in which the 
(spatially uniform) temperature T(t) of the JTJ varies in time as in (||). 

We reiterate that, in our approximation, the Zurek-Kibble bounds rely 
on the fact that the Swihart velocity is, numerically at least, the maximum 
speed at which the order parameter can change at any time. For speeds 
lower than this, we retain the adiabatic approximation, in which the corre- 
lation length of the field Xj(t) is just Xj(T(t)), where Xj(T) is the adiabatic 
(equilibrium) coherence length at temperature T. 

For a symmetric JTJ 



Aj(T) = V We(T)J c (T)' (14) 

where J C (T) is the critical Josephson current density introduced earlier and 
d e (T) is the magnetic thickness. As for the latter, if Xl(T) is the London 
penetration depth of the two (identical) superconducting sheets, 

then 

d e (T) = d ox + 2 X L (T) tanh - - %- . (16) 

2Al(J j 
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Neglecting the barrier thickness d ox <C d s , Xl, gives d e = d s close to T c , i.e 
the magnetic thickness equals the film thickness and can be set constant in 



All the T-dependence (and t-dependence) of A j resides in J c , which has 
the formlll 

jj^'m^m. fir, 

y ' 2 e PN 2k B T y ' 

In (|l7[) , A(T) in the superconducting gap energy and varies steeply near T c 
as 



A(T)~1.8A(0)Wl-^ (18) 

and pm is the JTJ normal resistance per unit area. 

Introducing the dimensionless quantity a = 1.6A(0)//sbT c , whose typi- 
cal value is between 3 and 5, enables us to write J C (T) as 

Jc(r)~aj c (o)(i-^). (19) 

J- c 

Thus, on repeated substitution of T(t) of @ in the above, we see that in 
the vicinity of the transition, 

\j(t) = \j(T(t)) = Z 0] ff (20) 

where 



V2 £/ M s aJ c (0)' (21) 

At the same time, for a finite electrode thickness JTJ, the Swihart velocity 
takes the form 



dox 



where 



di(T) = d ox + 2\ L {T) C oth(^^). (23) 
In the vicinity of the transition it follows, for the quench (0), that we recover 



c(t) of (PT[), c(t) = coJt/TQ, where cq = cqJ d s d ox /e r \ L {Q). Inspection of 



13|) shows that 



h = (24) 
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with tq = fo/co > in analogy with (^) . 

Making the second assumption of Zurek and Kibble that \j(tj) also 
characterises kink separation at this formation time, we find 

Aj = Xj(tj) = ivM = £o(^)*, (25) 

V t j Tq 

identical in form to (0) for superconductors. The quarter-power dependence 
of Aj on Co (for given £q) that we alluded to earlier is apparent. Moreover, if 
we keep c~o fixed, then Aj oc ^ 4 Inserting reasonable values of £o > 10 /xm 
and Co = 10 7 m/s gives f > lps. Thus, if tq =1s we have tj = 1 fis. We 
see that t s < tj by up to a factor of a few. 

Whatever, we have a practical problem with symmetric JTJs in this 
parameter range, in that substituting these values in Aj of (25) gives Aj « 
10 mm. This Aj which, in Zurek's picture, should characterise fluxon sepa- 
ration at the quench, would require very large annuli in order to see fluxons. 

Fortunately, such concerns are theoretical since the fabrication of JTJs 
typically yields non-symmetric devices with more acceptable properties. 



4. NON-SYMMETRIC JTJs 

Suppose the two superconductors, 1 and 2, now have different critical 
temperatures T C2 > T Cl . Fluxons only appear at temperatures T < T Cl , from 
which we measure our time t. At this time the individual superconducting 
gap energies are 



T, 



A 2 (T Cl )~1.8A 2 (0)Jl--^ 



(26) 



L ('2 



and Ai(T) vanishes at t = as 



A 1 (T C1 )^1.8A 1 (0)J— . 



(27) 



The critical Josephson current density J' c (T(t)) for a non-symmetric JTJ is 
nowlia 

7rAi(T)A 2 (r) ^ fr ,. A ,^! r a . .9^n-V2 



J'dT) 



(3{T)ep N 



E (K + A?(T)] [cf + Ai(T)]}- 



(28) 



/— — oo 



where /3 = 1/ksT and the = {21 + 1)tt/(3 are the fermionic Matsubara 
frequencies. Near t = (T = T C; i) 



J&T{t))~c[J'M\l 



T, 



t 



T C2 V T Q 



(29) 
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instead of (Eq), where 



cK ' [Ai(0)+A 2 (0)]ep^v' k B T Cl K ' 

provided A2(T Cl ) <C 2irkBT cl . This is the case here. 
If we now construct Xj(T(t)) we find 

Xj(t) = Xj(T(t)) « f (l - ^)~ 1/4 (^) 1/4 , (31) 



where £o is as in (pi), and we have used the fact that Jc(0) ~ J c (0) and 
a' ~ a. The crucial difference between A j(t) of ( |3l"| ) and Aj(t) of ( |20| ) for 
the symmetric case lies in the different critical index. \j(t) of ( |3T| ) is very 
insensitive to the time at which fluxons form. For the critical behaviour of 
( p9| ) to be valid, our earlier condition now becomes 1— T Cl /T c , 2 S> 0(Jj/tq) = 
O(10 -6 ), which is always the case. The critical time tj, determined from 
fl24| ) as before, now gives 

h = ffr%\l - (32) 

rather than (^J). For a typical value (1 — T Cl /T C2 ) = 0.02 and the same 
values of fo and tq as used previously, we find tj 0.24 /is. In turn, this 
gives 

Xj = Xj(tj) * &(1 " ^)~ 1/4 (^) 1/7 ^ 1.4mm. (33) 

This is an order of magnitude smaller than in the symmetric case, but still 
an order of magnitude larger than the corresponding lengths in annuli of 
single superconductors. 

However, we note that tj of (^) is now slightly smaller than t s of 
For all our caveats about the order of magnitude nature of the bounds, if we 
replace tj by t s , Xj of ( |33"|) is reduced slightly to 1.1mm. The difference is 
negligible, given the crudity of the bounds. 

In summary, for annular JTJs of this type we expect A j ~ 1 mm, rather 
than 10 mm or 100 /xm, and we would expect to start to see fluxons once the 
circumference C ~ 1 mm. This is, indeed, what happens. 



5. PAST AND FUTURE EXPERIMENTS 

As we said earlier, the original experiments^ performed by one of us 
(RM) were devised with other aims in mind. The intention was to produce 
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fluxons for further experiments and the density and frequency with which 
they were produced was secondary. The experiments that have been per- 
formed, and which we have in mind, are simple in principle. As we cool an 
annular JTJ uniformly from above its critical temperature T c , domains in <p 
will form whose boundaries are the fluxons. The mean flux is zero. 

The relevant observational quantity here is the variance in the flux or, 
equivalently, the variance AN in the number of fluxons minus antifluxons. 
This stability gives us a distinct advantage over experiments in vortex pro- 
duction in superfluid 4 He inwhich what is measured is a (decaying) density 
of vortices plus antivorticesEj'Ej. In fact the null result of Ref.|l2] for finding 
vortice s in 4 He can be attributed to uncertainty over the decay rate of vortex 
tangles^. 

If the circumference of the annulus is C = 2nr, and the mean tp corre- 
lation length at the time tj of the fluxon formation is Xj then, in parallel 
with our earlier observations, we would expect, for C S> Xj, 




of = ° UZ1 T- , (34) 



as a result of a random walk in phase around a ring. It is the prediction 
<M) that we would like to test for a variety of junctions and temperature 



quenches, with different Xj. Unfortunately, for the experiments to date it 
has not been possible to have C S> Xj. 



For the experiments of Ref. 1C non-symmetric annular Nb/Al-A10 x /Nb 
JTJs were fabricated with a whole-wafer process in which the junctions are 
formed in the window of the SiO insulating layer between the base and the 
counter electrodes. See Ref.|lO| for details. The ring-shaped junctions had a 
mean radius f = 80 /im, width Ar = 4 fim and a geometry in which both 
the base and top electrode have a hole concentric to the ring. The junction 
temperature was raised above its critical value (T c ~ 9.2 K) by means of a 
heating resistor placed close to the (unbiased) sample; then the temperature 
was lowered by letting the sample cool down by exchanging heat with the 
liquid helium bath through some helium gas. By changing the helium gas 
pressure it was possible to vary the thermal constant and hence the quench 
time tq in the range 10 -2 s-ls. The samples were measured in a very well 
electrically and magnetically shielded environment. The number of fluxons 
trapped during the transition was determined simply by feeding a proper 
current to the junction and measuring its voltage. 

Two different samples were produced, both with 1 — T Cl /T C2 ~ 0.02. 
For sample B, we estimated £0 = 6.5 ^m, cq = 10 7 m/s and tq = 0.65 ps. 
Although tj — t s , if we take ((p3|)) at face value, we find Xj ~ 1 mm (with 
experimental uncertainty of up to 50%). [ On the other hand, if we estimate 
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A j at Xj(t s ) we find Aj ~ 0.7mm, identical at our level of approximation.] 
With C ~ Aj we cannot use as it stands. However, we would 
expect to see a fluxon several percent of the time. In practice (invariably 
single) fluxons formed once every 10-20 times. 

In Sample A, having a similar circumference and similar superconduc- 
tors coupled differently, whereby Aj was 2.5 times larger, there was much less 
fluxon trapping, with no reliable fluxons seen in several quenches (enough 
to have seen something in the former case). This demonstrates that the 
superconductors cannot be treated individually for fluxon production, as 
our earlier discussion confirms. We are imprecise in the quantitative details 
since the aim was not to count the frequency with which defects occurred, 
but just to have defects at all. Nonetheless, it is apparent that, in order to 
have Aj <C C for annular JTJs of a size comparable to the one above, we 
must find ways to reduce Aj further, as well as possibly increase C. 

The quench time tq is difficult to establish accurately. Fortunately, 
the small power of tq/Tq in ( p3| ) makes the uncertainty largely irrelevant. 
However, it has the consequence that a huge reduction in the quench time is 
necessary to have any observable effect on Aj. Similarly, we can gain little 
from increasing the asymmetry of the JTJ, even if that is easily possible 
experimentally. In consequence, the only realistic way to reduce Aj is by 
reducing £o of (|H)- Primarily, we need to increase J c (0) to the largest 
possible value which does not degrade the barrier quality, J c ~ 10 4 ^4/cm 2 . 
However, for such currents we shall have i s > tj. The simplest assumption 
is that it will be i s that will set the time-scale at which fluxons are formed. 
New experiments with such values are under active consideration by us at 
this time. 

We consider this an important test of the Zurek-Kibble programme for 
establishing common constraints in the production of defects in condensed 
matter physics and relativistic quantum physics, given that the spontaneous 
creation of flux has not been seen in high IVplanar superconductors at the 
level suggested by the Zurek-Kibble boundEJ (in fact not at all), and that 
vortices arenot observed in 4 He quenches. This is even though a separate 
experimented supports the phase separation of the Kibble mechanism. 

We conclude with one final concern. Although the absence of fluxons in 
our Sample A, with larger Xj/C, suggests that the observed fluxons are not 
artefacts of the heating and cooling environment, we might be concerned 
at the effects of temperature inhomogeneities. In general, moving phase 
boundaries leave less (or no) defects in their wakeH. 

The critical slowing down of c{t) and c{t) weakens the effect of any large 
scale inhomogeneities in temperature, even though it is modified slightly for 
idealised non-symmetric JTJs. It is difficult to put quantitative limits on the 
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permissable inhomogeneities but, with empirically comparable £o 5 To, and tq, 
the situation is no better of worse for JTJs than for any other superconduct- 
ing system undergoing a mechanical quench. A possible check of these ideas 
could be performed by cooling with a temperature gradient across the ring 
circumference, although this would cause its own experimental problems. 
We postpone such considerations to the far future. 
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